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Abstract 

We study evolution of bipartite entangled quasi-Bell states in a qubit-oscillator 
system in the moderately strong coupling regime and also in the presence of a 
static bias. Using the adiabatic approximation that relies on the separation of 
the time scales of the two components, the reduced density matrices for the qubit 
and the oscillator degrees of freedom are obtained in closed forms comprising of 
linear combinations of Jacobi theta functions. The analytical results based on theta 
functions are found to be in good agreement with their series counterparts in the 
said parametric regime. The entropy of the system quantifying the entanglement 
is computed via the qubit reduced density matrix. The time dependence of this 
density matrix owes its origin to the bias term. The reduced density matrix of 
the oscillator is employed for obtaining the phase space distributions such as the 
Husimi Q-function, which, in turn, is utilized for extracting expectation values of 
antinormally ordered operators. 
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I Introduction 



One of the simplest and most important models in quantum physics describes a two- 
level system (qubit) interacting with a simple harmonic oscillator. The system is well- 
understood [T] in the context of the rotating wave approximation that holds for a small 
detuning i.e. near-equality of the oscillator and the qubit frequency, and a simultane- 
ously small qubit-oscillator coupling compared to the oscillator frequency. A property of 
this approximation is the conservation of the total number of excitations. In this regime 
the system is exactly solvable [1] allowing closed-form evaluation of dynamical quantities. 
Recently, however, a variety of experimental situations pertaining to stronger coupling 
domain, where the rotating wave approximation is no longer valid, have been investi- 
gated. Superconducting qubits interacting with nanomechanical resonators [2j|3] is one 
such promising area where experimentally the qubit-oscillator coupling constant may be 
driven to values comparable to the qubit and oscillator energy scales. Incorporating terms 
in the qubit-oscillator interaction Hamiltonian the authors of Ref. E] advanced an adi- 
abatic approximation scheme that utilizes the separation of the time scales of the qubit 
and the oscillator. The djTiamical state of the fast-moving component, say the oscillator, 
is assumed to instantaneously adjust to the slow-changing state of the qubit. This allows 
for decoupling of the full Hamiltonian into sectors related to each time scale, and approx- 
imate eigenstates of the system may be derived. The problem has also been studied by 
using a generalization of the rotating wave approximation [7] where the energy eigenstates 
have been investigated. For a single qubit analytical expression of the time- dependent 
behaviour of the two level system was obtained [5] in the small coupling limit. Studying 
a system of two qubits simultaneously coupled to a single oscillator mode the authors of 
Ref. [8j analytically evaluated the collapse and revival dynamics in the moderately strong 
coupling regime. Their procedure involved using Poisson sum formula [3] that leads to 
expressions of the reduced density matrix of the qubits as an infinite sum. 

It has been noted in [6l |T0] that whereas in the case of cavity electrodynamics the usual 
experimental set up is described by an unbiased qubit operated through its degeneracy 
point, the static bias of a superconducting qubit may be easily varied, say, by operating a 
magnetic flux on a Josephson junction. The energy spectrum of the system in the strong 
coupling regime in the presence of a bias has been investigated [TU] by employing the Van 
Vleck perturbation theory [TT]. Using a semiclassical approach the ground state in the 
presence of a bias has been investigated in [6]. Recently the effective Hamiltonian of two 
superconducting qubits coupled to a oscillator in the presence of a static bias has been 
explored 

In the current work we investigate the moderately strong qubit-oscillator coupling 
constant regime in the presence of a dominant oscillator frequency and a static bias of 
the qubit. In particular, starting with a quasi-Bell bipartite initial state we use the 
adiabatic approximation to find the time-evolution of the reduced density matrices of the 
qubit and the oscillator, respectively. Retaining up to the quadratic terms in the Laguerre 
polynomials in the phase factors the density matrix elements may be analytically expressed 
in closed forms via the Jacobi theta functions. Using this structure we explicitly obtain 
the entropy of the reduced qubit density matrix, which measures the time-variation of 
the entanglement of the bipartite state. To evaluate the analytical expressions we treat 
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the bias perturbatively, and retain terms of second order in smallness. The summation 
procedure may, however, be extended to arbitrary order. Moreover, in the course of our 
analytical derivation we neglect the bias terms in the exponents of the phase factors as 
that will be of importance in the asymptotically long time limit. Here we restrict our 
attention to a finite number of revivals. The reduced density matrix of the oscillator 
yields the time evolution of the Husimi Q-function, and this phase space distribution is 
subsequently used to analytically evaluate averages of the antinormally ordered operators. 



II The reduced density matrices 

The boson-qubit Hamiltonian considered in our work reads, in natural units (h = 1), as 
follows: ^ 

H = ——ax — -o"2 + ua) a + Acr^ {a) + a), (2-1) 

where and a are the raising and lowering operators of the harmonic oscillator with 
frequency w, the Pauli spin operators and describe the qubit with energy splitting 
A as well as static bias e, and the qubit-oscillator coupling strength is denoted by A. 
The standard Fock states {|n), n = 0, 1, . . . ; a\n) = -^/n \n — 1), |n) = ^/n + 1 \n + 1)} 
provide the basis for the oscillator, and the following eigenstates cr^l + 1) = + | + 1) 
span the space of the qubit. The Hamiltonian fl2.1l) is not known to be exactly solvable. 
It departs from the rotating wave approximation that holds good in the near-resonance 
and weak coupling regime. In the present work we follow, as mentioned earlier, the 
adiabatic approximation [5J that relies on the separation of the time scales characterizing 
the high oscillator frequency and the (renormalized) low qubit frequency. The fast-moving 
oscillator then adiabatically adjusts to the slow changes of the state of the qubit. 

To proceed one notes that, in the said time scale, the oscillator follows the instanta- 
neous state of the qubit observable a^, and for the oscillator degree of freedom the spin 
variable az may be assumed to be in its eigenstate with the eigenvalue given by +1. The 
effective Hamiltonian of the harmonic oscillator now reads 

Ho = oja)a±\{a) + a). (2.2) 

The number states shifted symmetrically by the displacement operator diagonalize the 
above Hamiltonian: 

|n+> = (iA/w) |n>, T){a) = exp(Q;a^ - a*a), (2.3) 

where the degenerate eigenenergies are as follows: 

A^ 

E^^£n = nuj . (2.4) 

u 

The number states displaced in the same direction maintain orthonormality, whereas 
the overlap between the Fock states displaced in opposite directions is given by 



(—1)™ " exp (— I) x^z^^nl/m! " {x) , m^n 
(m_|ri+)= \ — — — (2.5) 
exp (— |j X 2 -^m\/n\ [x) , m < n, 
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where the associated Laguerre polynomial reads Li^{x) = Xifc^o (~-'-)'^ (n-fe) IJ' 
parameter x = (2A/a;)^ acts as the perturbative expansion coefficient. 

After diagonalizing the high frequency oscillator component of the Hamiltonian one 
now turns attention to the low frequency qubit part. Tensoring the qubit states with the 
displaced oscillator basis {| + l,r;,+)|VT2, = 0,1,...} the Hamiltonian may now be readily 
constructed. In the parameter regime A < a; one may neglect the matrix elements that 
mixes oscillator states with different eigenvalues of its number operator. This reduces 
the Hamiltonian to a block-diagonal form where each block mixes displaced oscillator 
states with identical number of photonic excitations. The Hamiltonian for the n-th block 
assumes the form 



H 



nu 



V 



nu 



he 



A ~ fx 

bn = -— L„(x), A = Aexp 



e 
2 



(2.6) 



that may be diagonalized in the basis of the eigenstates listed below 




where the corresponding eigenvalues read 



l,n_>, 



e2. 



(2.7) 



(2. 



In order to understand the qubit-oscillator entanglement we study the dynamical evo- 
lution of the quasi-Bell initial states of the coupled system in the presence of a bias 
(e^O): 

|^(0)>(±) = (|1, a> + I ^ 1, -a» , |a> = D(a)|0>, (2.9) 



where the coherent state {\a) Va = 3?(q;) -t- i$j(a) e C} for the oscillator degree of 
freedom may be realized by the action of the displacement operator on the vacuum. 
The nonclassical entangled states (12. 9 p have been discussed earlier [13] in the context of 
generating bipartite and multipartite entanglements in cavity quantum electrodynamics. 
It has also been observed [6] that in the large coupling regime (A > oj) a state of the 
generic type (12. 9 p becomes the approximate ground state of the combined system. In 
the present work we, as an example of a typical nonclassical state, analytically study 
the properties such as time evolution of its entanglement measure and the phase space 
distributions of the oscillator for moderately strong qubit-oscillator coupling. The precise 
limits of the coupling constant and the bias parameter will be made explicit later. The 
time evolution of the initial state (12. 9p is given by 



\m) 



(±) 



V2 



exp 



00 -. 
n=0 



a 



H±exp(-2Xn^)IO + 



, (2.10) 
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where a = a + X/u, = and = a/^^ + (~1)"tT"Ta/ ^ ■ ^'^^ 



If^ral V '^X.r, 

f l2.10p and hereafter the upper and lower signs refer to respective quasi-Bell states given 
in (12. 9p . Employing the state (12.101) the evolution of the density matrix of the composite 
bipartite system may be readily constructed: 

p(±)(t) = \m)^-\m\- (2.11) 

The marginalized density matrices for the individual degrees of freedom may be obtained 
by partial tracing of the composite density matrix (12. lip . For instance the tracing of the 
oscillator degree of freedom yields the reduced density matrix of the qubit. Its explicit 
structure obtained via (I2.11[ I2.10p assumes the form 



pS|Hi)-Trop(±)(i) 
where the components read 



(2.12) 



00 2n • 2 J. 

C - e~exp(-,l)2(-ir^5.^^, (2.13) 

n=0 

e = - exp(-^s)2j(-l) l-2e +^^^ <^-l^+>- (2.14) 

^ „=0 ^ ^'^ An / 

An interesting feature evident in the reduced density matrix (I2.12p that it is time- 
independent in the limit e = 0. Physically one may conclude that when the fast-changing 
oscillator degree of freedom is traced out the slow-moving qubit degree of freedom is close 
to its maximally mixed form that becomes time-independent in the bias- free limit. 

Similarly the trace over the qubit basis generates the reduced density matrix of the 
oscillSiijOr ' 

p'~o\t)^^^Q\mf-\m\ (2.15) 

that may be expressed, via (I2.10p . in terms of the displaced number states as follows: 

2 „tlo vn!m! V 

+ (-l)"+™<+)*(t)C:W(t)|n_Xm_|) exp ( - ^(n ~ m)wt), (2.16) 
where the coefficients read 

= exp(zx„t) + exp(-zxnt), (2.17) 
^ X^ + i±(zl)!^^ ^,^,^ Xn-e±{-m^ ^ ^2.18) 

2Xn 2^n 

We note that apart from studying in the context of adiabatic approximation no other 
simplification has been made in the derivation leading to the reduced density matrix 
elements fl27[3l ICTl [216|l . 
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Ill The qubit density matrix 

We now focus our attention towards analytical evaluation of the density matrix in the 
parameter regime where the oscillator frequency sets the dominant scale and the following 
hierarchy is maintained: e < A < a;. In the weak coupling limit A « a; the series in fl2.13[ 
I2.14P may be summed by retaining d /a [S] only linear terms in the Laguerre polynomials: 
Ln{x) = l — nx + 0{x'^). For simplicity we retain terms O(e^) in the bias parameter. Since 
we are interested in the first few revivals of the system we also systematically neglect e 
dependent terms in phases as they will be important in very large time limit. The higher 
order terms in the bias parameter e may, however, be taken into account in the present 
perturbative scheme. In this parametric domain the elementary trigonometric identities 
now allow us to compute the diagonal and the off-diagonal elements of the density matrix 
given in (12.12112.1^ in the following form: 

/cosl^f)), (3.1) 

2f?|cos2^) X 

(3.2) 

where the scaled time variable as well as the phases are given, respectively, by r = 
At,^lt = q\ sin(xT) + r, Q!f^ = q\ sin(xr) + (1 — x) r, and the parameter / = xq\ deter- 
mines the amplitude of the mode with the renormalized frequency. This approximation 
that retains up to linear terms in Laguerre polynomials hold in the weak coupling range. 
For higher qubit-oscillator coupling the linear approximation method fails for the higher 
revival times. 

In an attempt for analytically evaluating the density matrix elements in the context 
of Tavis-Cummings model for higher values of the qubit-oscillator coupling constant the 
authors of Ref. ^8j| considered a method based on Poisson sum formula. Here one retains 
contributions of the order O(x^) ^ 0{{\/uY) that amounts to considering up to the 
quadratic terms in the expansion of the Laguerre polynomials: L„(x) = l — nx + ^ x"^ + 
O(x^). Following [8] we now provide a Poisson sum for the elements of the density matrix 
fl2.12p . To proceed, for a large value of we approximate the Poisson distribution in 
fl2.13[ I2.14P with the corresponding Gaussian limit: 

P(„).expH|)f .-±^exp(-i!i^). ,3.3) 

The Poisson summation formula employs Fourier series and, in the present context, it 
may be stated [9] as follows: An arbitrary function g(x) integrable in the domain (— oo, oo) 
satisfies the sum rule 

t ; - i: f%(t)exp(-2.«t)*. (3,4) 

n= — 00 £= — CO °° 

For the case » 1 the lower limit of integration corresponding to (13.31) may be extended 
to (— oo), and the integrals in the Fourier transforms may be evaluated by the standard 
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c 



^2A ~ '^^ ("2^1) + ^exp (^-2^1 cos^ yj (^cos^ 



2A 

1 + / 



cos Qt " f cos Q 



(x) 



2AA 

+ exp (^-2qI cos^ y ) sin Qt, 



Gaussian integration technique. With the above structure in mind we obtain the Poisson 
sum for the diagonal term (12.131) in the density matrix: 



c 



e 

2A 

e 



exp (-2^,1) (^1 - /(I + x) + ^-f^ 



£=0 



(3.5) 



In (13. 5 p the terms comprising the Fourier components read 
= exp ( — (j^-?) ( r — vl'^' ' 

J6{1,2,3} V 

i ^ <j I tan 



x; 



0) 



X exp 



l)+2v.J^)(r) 



-0) 



,(3.6) 



X f T 1 

where the parameters are given by tan^ = ^1^6% The strength of the 

Gaussian packets in the rhs of (13. 5p on the respective powers of the parameter /. The 
values of the scaled time variable r at which the peak of the contributions (13. 6p occur, 
and the decay of the corresponding exponentials are given, respectively, by 



V 



0) 



27r(£+ 1/2) 



-0) 



X 1 + 



(3 - 2j)x f 
4 2 

The phase constructs in (13. 6p read 



2 , (3 - 2j)x f 



(3.7) 



^?{r) = x{l+^^^-^)r-2n{i+l/2) 



{3.t 



and the scaled time-dependent parameters are given by 



r 



(1) 



r, r 



(2) 



;i - x)t, t 



(3) 



l-2x + 



X 



T. 



(3.9) 



We mention that in evaluating the first term in the rhs of (13. 5p the Gaussian limit (13. 3p 
has not been exploited as the corresponding term in (I2.13P is free from a complex Laguerre 
polynomial dependent phase and is amenable to elementary summation technique. Simi- 
larly, in summing other series expressions we will continue to use elementary procedures 
wherever applicable. 

Continuing the above Poisson sum process the off-diagonal element (I2.14p for the 
reduced density matrix of the qubit may be evaluated as 



e = ^exp(-2^|~x/2) fl + / + ^) - 



£=0 



2AA 
3 
4 



exp(-2£.|) (l-/(l+x) + ^/^ 



-v^ S ^[^i'^ - /(I + ' + -J'^e'^)] - ^ S ^i^i'^ (3-10) 



The contribution of x]"''' in the reduced density matrix elements fl3.5[ I3.10p . as evident 
from f l3.6p . is maximum when its real part achieves maximum value at r = v'f'^ that 
corresponds to successive revivals of the two-level system. The index therefore, has 
the natural significance of enumerating specific revivals of the system. A sharp peak 
corresponds to higher value of the fall-off parameter <^^-'''. From the structure of (13. 6p it 
follows that for a positive r only the Fourier components with positive indices £ ^ need 
to be summed over, as it was maintained in obtaining the density matrix elements (13. 5[ 
[MOD. 

In the said limit of moderately strong coupling regime where we retain contributions 
O(x^) in the density matrix elements and consider up to quadratic terms in the Laguerre 
polynomials an alternate method of summing the series in (12. 13^ I2.14p may be proposed 
by employing the Jacobi theta functions. This method has the additional benefit in that 
it expresses the matrix elements in closed form. We quote [9] the series expansions of the 
Jacobi theta function that is relevant in our case: 



d^{q,z)= g"'exp(2mz), ^M.z) = (-l)V'exp(2m2). (3.11) 



These interrelated functions d^{q^ z)^di{q, z) depend on two complex arguments and the 
series sums in (13. lip are convergent for the choice \q\ < 1. Our starting point for evaluating 
the infinite sums in the series (12.1311 2.14p via the theta functions is that we again use the 
Gaussian limit expressed in (13. 3p . For a large value of q\ the lower limit of the series in 
(I2.13[ 2.14P may be extended to (— oo). Due to the sharp decay of the Gaussian function 
the error encountered is quite small. This process is similar to the evaluation of Fourier 
components via the Gaussian integration technique in the previously described Poisson- 
sum method, where the lower limit of integration is extended to (— Qo). We here mention 
that a truncated series on one side in the definition (13. lip remains convergent in the 
domain |g| < 1, and may be viewed as a close analogue of the theta function. Since the 
difference between the alternate cases in our parametric regime is quite small we retain 
the standard definition (13.110 . The resultant closed-form expression for the diagonal term 
(I2.13P now assumes the form: 



o , exp 

exp( -2f,|)(l-/(l + x) + -f) 



5R[0(r)] , (3.12) 



where 0(t) is explicitly given as a sum of theta functions as follows: 

0(t) = exp (-2x^)^94(9,2^^^) "/(I + 2:) exp (-2r(2))^?4(^,z(2)) 

+ ^/2exp (-zr(3))^4(g,;2(3)). 



(3.13) 



The three theta functions appearing in the combination (I3.13P posses same value of the 
complex argument q but differing values of the other complex variable given, respectively, 
as 



g(T) = exp 



2^1 



exp I — I 



z^'Ht] 



Jg{1,2,3} 



l(.(l + (3-2,)9r 



(3.14) 
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For the off-diagonal element (I2.14p a similar closed-form expression involving theta func- 
tions may also be obtained as follows: 



x.l-/(l + x) + ^f)-'"P 



2AA 

m 



exp ( - 2^1) 



5J[e(r)] 



6 exp ( - g|/2) 
2A 



[exp(-2rW)i94(?,^^'^)]. 



(3.15) 



This completes our procedure of evaluations of the elements of the reduced density matrix 
of the qubit via the Jacobi theta functions. It involves considering up to the quadratic 
terms in the expression of Laguerre polynomials appearing in the arguments of the phase 
factors and also approximating the Poisson term with its Gaussian limit. For moderately 
large coupling (A/w ~ 0.25) this summation technique has good agreement with the 
numerical evaluations based on fl2.13[ 12.141) . Concerning our choice of bias parameter 
we note that as we have included terms up to O(e^), we may maintain the order of 

the dimensionless ratio — ~ 0.1. These restrictions describe the parametric regime in 

A 

which our analytical evaluations based on theta functions hold. In Figs. [T]|3] we now 
compare the closed-form expressions of the density matrix elements fl3.12[ I3.15P with 
the corresponding sums given in (I2.13^ I2.14p . For a significant number of revivals good 
agreement is observed. 

(a) 





(b) 





Figure 1: The diagonal component C of the density matrix for (i) A = 0.1a;, A = 0.15^, e = 
O.Olw and a = 2 (left) and (ii) A = 0.18a;, A = 0.3a;, e = 0.02a; and a = 1 (right) using 
(a) the series fl2.13p . and (b) the evaluation fl3.12p using theta function. 
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(a) 




X/27I; x/27r 

Figure 2: The real part of the off-diagonal component ^ for (i) A = O.lcj, A = 0.15a;, e = 
O.Olu and a = 2 (left) and (ii)A = 0.18a;, A = 0.3a;, e = 0.02a; and a = 1 (right) using (a) 
the series fl2.14p . and (b) the evaluation fl3.15p using theta function. 



The time evolution of the entanglement of the initially chosen quasi-Bell state (12.91) is 
of interest. For the bipartite system considered here this may be measured by the entropy 
of the reduced density matrix of the qubit: S{pq) = — Tr(pg Inpg). The construction 
(I2.12P of the reduced density matrix pq produces its pair of eigenvalues given below: 

A± = ^ + VC^TW, (3.16) 



which are valid for both the initial choices (12. 9 p of quasi-Bell states. The said entropy 
now reads 

5 = -A+logA+ - A^logA". (3.17) 

The structure of the diagonal and the off-diagonal terms given in (I2.13[ I2.14p may now be 
employed to extract the eigenvalues of the density matrix (12.121) up to the order O(e^): 

A± = i (l + exp{~gl - x/2) Ve) , (3.18) 

where the sum is given by 

00 2(n+m) 

6 = V {-l)(-+"^)^^^L^{x)LUx) X (3.19) 

7i,m=0 

r ^ sin^iXnt) 2 ^M'^Xnt) sm{2xmt) ^/^Ax 2 sin^(xn0 sin^(xmO ] 
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(a) 




Figure 3: The imaginary part of the off-diagonal component ^ for (i) A = 0.1ci;,A = 
0.15u;,e = O.Olw and a = 2 (left) and (ii)A = 0.18w,A = 0.3a;, e = 0.02a; and a = 1 
(right) using (a) the series fl2.14p . and (b) the evaluation fl3.15p using theta function. 



It is interesting to consider the entropy in the bias- free e = limit. Towards this end we 
use the identity 



y ^^L„(x)t" = exp(-x)/o(2V^), (3.20) 
where the Bessel function is denoted as 

(-1) 

The equations ranging f l3.17tl3l2T|l may now be used to obtain an exact closed-form ex- 
pression for the entropy in the bias- free limit: 

S |,.o = ln2 - ^ exp ( - 4^1 - x) {l,{2^f))\ (3.22) 

The entropy fl3.22p in the e = limit is independent of time and attains its maximal value 
for a large Qa- 

It is evident from f l3.19p that the inclusion of a bias term modifies the entropy (13.221) 
only at the order O(e^). We first obtain the entropy in the weak coupling limit where the 
linear approximation of the Laguerre functions as advocated in [S] is assumed to be valid. 
Our previous constructions given in (13. 1[ 13. 2[ I3.17[ I3.19P now produce the entropy in the 



10 



said limit: 



S = ln2-iexp(-4f5|-x)[l + 2/-|^[l-((l + /)cosfi,-/cosfiS^)) x 
X exp (2gl sin^ - ^ exp (igl sin^ y ) sin^ ^* ~ ^ exp(x)(l - 2/) 
-^exp(^4^|sin2 Y + x)(^cos^f]t-2/cosacosfif^) j j . (3.23) 

It is easily evident that the expression of entropy (13.231) in the weak coupling limit agrees 
in the bias-free e = case with the evaluation (13.221) when the rhs of the latter is expanded 
retaining up to terms linear in the parameter /. 

For the moderately strong coupling described by the parametric regime specified in 
the context of the equation (I3.15P we evaluate the entropy (13.171) by retaining terms up 
to the order O(x^) in the Laguerre polynomials. The summation procedure as evidenced 
in (13.121 13.15P now yields the entropy via linear combinations of Jacobi theta functions. 
Up to the order in the bias parameter O(e^) the entropy reads 



S 



ln2 - - exp {^Ag, 



X) 



l + 2/ + ^f -S^[l-/x-A-(l + / + 



4 



exp (3^1/2) 



+ 



exp {3gl/2 



5R[e(r)] 



exp(3gg 



?i/2 + x) 



[exp {-iT^^^)di{q,z^^^ 



exp [?,g\ + x) 



5R[e(T)]) - 



exp(x) 



1 - 2/(1 + a;) + r - + 2x 



(3.24) 



We present below in Fig. [3] the plots of the time evolution of the entropy of the system 
based on our closed- form expression (13.241) and compare it with sum given in (I3.17[ I3.19P . 
For first few revivals good agreement is obtained. As physically expected the entropy dips 
from its maximal value during successive revivals of the system. The extent of the dip 
increases with the increasing bias parameter. For the bias-free case (e = 0) the entropy 
remains constant. 



IV The Q-function of the oscillator density matrix 

The Husimi Q-function [14j is one of the distribution functions that facilitates charac- 
terization of the quantum state of the oscillator degree of freedom. For our choice of 
quasi-Bell initial states (12. 9p it is defined as the expectation value of the reduced density 
matrix of the oscillator (I2.16P in an arbitrary coherent state: 

Q^^\P) = Up\p^o\tm- (4.1) 
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(a) 



CO 0.6915 




(b) 



CO 0.6915 




CO 0.689 




CO 0.689 



0.68V t 




Figure 4: Evaluation of the entropy for (i) A = 0.1a;, A = 0.15a;, e = O.Ola; and a = 1 
(left) and (ii) A = 0.18a;, A = 0.3a;, e = 0.02a; and a = 1 (right) using (a) the series fl3.17p . 
and (b) the evaluation fl3.24p using theta function. 



Substituting the said reduced density matrix evaluated in fl2.16p the Q-function of the 
state under consideration may be written as, 



^-p(-£,|) 2 (exp(-4) <°^'';,iT^''" (e!."'(i)e!,?'(t) 

n,m=0 



n\m\ 



(S/3*)" (S*^)^ 



+ (-l)"+™exp(-f;|; — 

X exp(— i(n — m)ujt), 



(4.2) 



where we abbreviate: (3 = (3 + X/tu, = \f3\ and ^ = j3 — A/o;, = \^\. To evaluate the 
sum in (14. 2 p we adopt a procedure discussed in the context of (13.11 13. 2p that, following 
[5], entails retaining up to linear terms in the expansion of Laguerre polynomials. As the 
two sums in (14. 2 p factorizes they may be easily obtained in the said limit. The positive 
definite Q-function now assumes the form 



= ^exp(-^>l)(exp(-4)|X(±)r + exp(-4)|2) 
where the complex density functions read 

X(±) = exp($l) cos ($,^ - ± ^ exp(-$!) sin (<l>,^ + ^ 



(±)|2 



(4.3) 



(4- exp($|) ^ sin ^$ 



i-]^+x^t sin (^^i 



)] 



+^- 



2A2 



exp(-$f ) \^ sin ^<l>^ ^ ^ ^ 



+ 



], (4.4) 
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(a) 







Figure 5: The Q''^\f3) function using (i) the series (14. 2p . and (ii) the analytical result 
(14. 3p . in left and right columns respectively, for A = 0.08a;,A = 0.15ijj,e = O.Olu and 
a = 2 at various time t (a)0 (b)300 (c)500 (d)900. 



and 

2)(±) = exp(-^!) cos + ^exp(^!) sin (^^^ - ^ 

-tj exp(-^!) [ sin + _ x sin + ^^—^) ] 

exp(vl>j) [ sin (v^,^ - ^) + sin (vp,^ - ^^^) ] • (4-5) 

In (14.41 145|) we have used the notations $t = 3^* exp(-ia;t), \E't = 3^* exp(-zcut), = 

XT II XT XT XT 

COS — , = COS — , $^ = $j sin — and \E'^ = \E't sin — . In the parametric region 

X/u < 0.1 where the linear approximation holds we now compare the analytical evaluation 
(14. 3 p of the Q-function with the corresponding double sum (14. 2p . The time evolution of 
the Q-function is faithfully reproduced by the closed- form expression (14. 3p . 

One of the utilities of the phase space distribution Q is that it provides a convenient 
evaluation of the expectation values of the operators expressed in their antinormal ordered 
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form. For instance, the first and second moments of the scaled electric field operator 
defined as 

^ (a-a^) (4.6) 



'2i 

may be expressed [H] via the following phase space integrals over the complex plane: 

<Ey±) ^ Tr(E(a,at)p(,±)(t)) = ^ (/? - /?* )Q(±)(/?)ci2/3, (4.7) 

<E2>(±) ^ Tr(E2(a,at)p(,±)(t) = ~-J ^((/3-/3*)2 + l)g(±)(/3M'/?- (4-8) 

Employing the Fourier sum fl4.2p and using the following notation for the coefficient 

= + {-^Y^^SmfKt)) exp(-^£a;t) (4.9) 

the first two moments of the electric field is observed as 

{^r-^ - 4^exp(-,i)2^5>[^Swl, (4.10) 

n=0 ^ 

<E^y±^ - el + \{i~^M~QV)Y.'-^^[G^kt)]). (4.11) 

?i=0 

With the purpose of evaluating the moments (14.101 14.111) in a closed form we now quote 
explicit values of the first two Fourier components given in (14. 9p : 

Sn;l{t) = ae(+(-l) {6n-6n+l) 

^ Xn Xn+1 

~i cos x„+it - cos exp(-za;f), (4.12) 

^ Xn Xn+1 J ' 

y;.2(t) = 2a I COS X„t COS Xn+2i + (e + dndn+2) 

^ Xn Xn+2 

+z(-l) [dn cosXn+2t- On+2 COS x„t exp(-22a;t). (4.13) 

^ Xn Xn+2 

Utilizing the definitions (13. lip and the construction given in (14. 8114. 13]) the first two mo- 
ments of the electric field as linear combinations of of Jacobi functions. The first moment 
is non- vanishing only in the presence of a bias term (e 7^ 0) in the qubit Hamiltonian (12. ip : 

<E>(±) = _4 e^P(-gi/2 ) /c^ (- ,aexp(-2a;t)] $^(T] + 55 (S exp(-za;t)] 3ft(H]y (4.14) 



A V47r^a 

where the linear combinations of "i^s and "di, functions read, respectively, as 

T ^ ((2 + x + x^)+(2 + |i)/exp(^)+^/^exp(^)) X 

X exp (^^)^93(q,z) + x(l + x) exp - %r(\ - |))t?3(g, z^^^) 

+ ^x/exp ( - ^r(\ - f + ^)) ^3(g,z(^)), (4.15) 
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and 



H = + x) - ^x/exp - exp (^^)T?4(q,z) 

+ x) exp - iT(^l - |j j?94(g,z(^)) 

+ ^a:/exp ( - zr{l - ^ + ^))Uq,z^'^). (4.16) 

The complex variables associated with the 'd functions listed in ( 14.15^ I4.16P are given 
below: 



Analogous to the above summation procedure based on d function a closed-form ex- 
pression may also be obtained for the second moment of the electric field (14. lip : 

+^ [la^ exp{^2tut)]'^[n]^ , (4.18) 
where the respective linear sums of and i}^ functions are abbreviated as follows: 



A 



exp (^iTx(^l ~ ^))^93(q,3) 



7^ 



+ ^exp (- tr[l - X + ^jjUq,i^)), (4.19) 
[ (1 - -^^(l + 2x + i^)) exp {^rx{l - ^)) + ^/(l + Ax) x 



2A2 V 2 // V V 4// A2" 

xexp (^irx(^l-^)) --^fexp (^zrx(l-^)) j i?4(q,3) 

+ ^x/exp ( - zr(l - 2x + ^))^4(q,3(2)). (4.20) 
The complex variables of the 1!} functions quoted in (I4.19[ I4.20p are listed below: 

This completes our demonstration of evaluations of the expectation values of antinormally 
ordered operators in closed forms comprising of the linear combination of i3 functions. 
The procedure is valid for the moderately strong coupling regime where we retain up to 
quadratic terms in the Laguerre polynomials. In Figs. [B] and [7] we compare, respec- 
tively, the expectation values (14.141 I4.18P with the corresponding sums (I4.10[ 14. lip in 
the parametric domain referred to in the context of (13.150 . We observe good agreement 
between two sets of values. 
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(a) 



(b) 





Figure 6: First moment of the electric field operator for the values A = 0.16 a; and a = 1 
derived by using (a) the series fl4.10p . and (b) the analytical evaluation f l4.14p . 



(a) 



(b) 



20 
T/27r 



20 

X/27I 



Figure 7: Second moment of the electric field operator for the values A = 0.16 a; and a = 1 
derived using (a) the series (14. lip , and (b) the evaluation fl4.18p using theta function. 

V Conclusion 

We investigated a qubit-oscillator bipartite system in the framework of adiabatic approx- 
imation in the regime of moderately strong coupling and also in the presence of a static 
bias term in the Hamiltonian. Starting with an entangled quasi-Bell state we obtained 
the evolution of the reduced density matrices of the qubit and the oscillator. The re- 
duced density matrix of the qubit is time independent in the absence of the bias term 
in the Hamiltonian. Identifying the series with the Jacobi theta functions we evaluated 
the density matrix elements in closed form as linear combinations of theta functions. The 
reduced density matrix of the qubit, in turn, yields the entropy of the system. On the 
other hand the reduced density matrix of the oscillator provides a way to study the phase 
space density functions of the oscillator degree of freedom. In particular, we evaluate the 
Q-function and utilize this for obtaining closed-form expressions of the first two moments 
of the scaled electric field operator as linear combinations of theta functions. The first 
moment of the electric field in the quasi-Bell states assumes non-zero value only in the 
presence of the bias parameter. Our closed-form evaluations of various physical quantities 
are compared with, and found to be good approximations of, their series values in the 
regime characterized by moderately strong coupling as well as weak bias . 

The theta function based summation procedure advanced here may be of utility in 
developing the theory for still larger qubit-oscillator coupling regime. In these instances 
we need to include higher order terms in the expansion of Laguerre polynomials appearing 
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in the phase factors of the density matrix elements. It is evident from our derivation that 
the corresponding series still remains convergent owing to the Gaussian weight factor. 
A generalization of theta functions depending on more variables may be developed to 
evaluate such quantities in closed form. Another natural direction in which the present 
approach may be extended is to study a system with a higher number of qubits coupled 
to an oscillator. These works will be presented elsewhere. 
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